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Abstract 

Measured response functions and low photon yield spectra of silicon photomultipliers 
(SiPM) were compared to multi-photoelectron pulse-height distributions generated 
by a Monte Carlo model. Characteristic parameters for SiPM were derived. The 
devices were irradiated with 14MeV electrons at the Mainz microtron MAMI. It 
is shown that the first noticeable damage consists of an increase in the rate of 
dark pulses and the loss of uniformity in the pixel gains. Higher radiation doses 
reduced also the photon detection efficiency. The results are especially relevant for 
applications of SiPM in fibre detectors at high luminosity experiments. 
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1 Introduction 



After a few years of R&D, silicon photomultipliers (SiPM) are today an in- 
teresting alternative to conventional vacuum phototubes in many applications 
due to the specific advantages of a solid state device: small size, low volt- 
age operation and magnetic field insensitivity 

a a a 

. Nevertheless, the idea 

of grouping hundreds of miniature avalanche photodiodes (APD) in a planar 
array to form an analog device has found a limitation in the high rate of in- 
trinsic dark pulses of such a configuration. Only small active area devices are 
currently available and typical noise rates are of the order MHz/mm 2 with 
signal pile-up due to optical cross-talk. This leads to the situation of a small 
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photo-sensor that needs to detect a relatively large amount of light in order to 
achieve an acceptable signal-to-noise ratio. To our knowledge there has been 
no attempt to use SiPM in applications where only a few photons must be 
detected. This is for example the situation faced when using scintillating fibres 
as tracking detectors. 

There is little information about radiation hardness of SiPM. It is well known 
from other types of silicon detectors that generating centres are created during 
irradiation which increase the leakage current [4{] . The bulk leakage current is 
multiplied in APD by the gain factor and the resulting pulses are undistin- 
guishable from photon generated events. Consequently, an increasing rate of 
dark pulses as a function of the radiation dose is expected for SiPM. Low light 
level detection will be degraded or even impossible in some applications if this 
effect happens to be large. Adverse effects of irradiation on other characteristic 
parameters of SiPM such as gain uniformity, after-pulsing or optical cross-talk 
probability will be also detrimental for a detector. In general, it is mandatory 
to study the impact of the particular kind of radiation the detector will be 
exposed to on the characteristic parameters that must remain stable. 

In this paper a complete model for multi-photoelectron pulse-height distribu- 
tions (MPHD) was used to extract the main consequences of electromagnetic 
irradiation on SiPM. Section 2 discusses the MPHD model. Section 3 gives 
a brief description of the irradiation set-up. Section 4 describes the changes 
in SiPM characteristics obtained for low and high irradiation doses. Section 5 
closes with future prospects of using this kind of detectors in low light level 
applications. 



2 The Monte Carlo model for the detector output 

A complete model for the MPHD of SiPM was derived based on a Monte Carlo 
simulation. For the model, Poisson distributed photo-electron and dark signals 
are generated independently. Each of these signals can cause optical cross-talk 
and after-pulses according to a given probability distribution. For this set of 
pixels the distribution in time is generated with respect to an integration time 
window. Single pixels can contribute with varying gain to the generation of the 
detector output. Finally, the noise is added. The free parameters of the simu- 
lation are: mean signal amplitude {A) ) gain variation (og), dark pulse rate (r), 
optical cross-talk probability (p pt), after-pulse probability (p a ft), trap lifetime 
(r), mean number of detected photons (A), pedestal position (x pe d), and noise 
amplitude (a pe d)- In the following, the simulated processes are explained. 

As SiPM are a set of APD connected in parallel many of their properties 
are thus inherited. Its dark count rate, r, is the sum of all the APD dark 
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Fig. 1. The simultaneous signals from up to three pixels are visible in this oscillo- 
scope trace. Optical cross-talk between adjacent pixels that form a cluster generates 
these events. The MPHD model described in the text assumes all signals in a cluster 
to appear at the same time. 

count rates. After-pulses appear in each pixel after a photon (or a thermically 
generated charge carrier) triggers an avalanche due to trapped carriers that 
are released after some time-delay. The after-pulse probability, p a ft, is defined 
as the fraction of events in which one additional signal is generated after the 
detection of a photon. It is the product of the trapping probability and the 
triggering probability, both increasing linearly with the applied bias voltage. 
The trap lifetime, r, determines the typical time-scale. 

SiPM are manufactured so that signal uniformity from pixel to pixel is quite 
good, typically within 10 % [5]. The small variation, ac, together with the nar- 
row single electron response function of each APD provides excellent photon 
counting capabilities and as many as 20 photons can easily be resolved in a 
typical pulse height spectrum [1]. The distance between multi-pixel peaks in 
the spectrum is a measure of the charge gain, A. 

Surface leakage currents do not cross the multiplication region. Their fluctu- 
ations are merged with other noise sources to define the noise amplitude as 
measured by the pedestal width, a pe d. A change in the surface leakage current 
should also appear as a shift in pedestal position, x pe d. 

An avalanche of 10 6 carriers in any of the micro-metric APD forming the 
SiPM will create around 50 photons via hot carrier luminescence with enough 
energy to trigger any neighbouring pixel [6]. Devices without trenches filled 
with opaque material exhibit optical cross-talk, where at least one photon 
is able to cross the spacing between micro-cells to produce a simultaneous 
signal (within 100 ps). Fig. [1] shows the effect of optical cross-talk in the SiPM 
elementary signal. The well defined amplitudes of the multi-pixel events is a 
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Fig. 2. (a) Measured single-electron response function of a SiPM. The multiple peak 
structure is to a large extend due to the high rate of dark pulses. Optical cross-talk 
and after-pulses contribute to a lesser extend, (b) Integration of noise signals in a 
SiPM with a randomly generated gate. Events outside the pedestal peak are a result 
of the high rate of dark pulses. 

result of the simultaneity of the composing signals. 

To the knowledge of the authors there has been only one serious attempt 
to describe the statistics of multi-photoelectron pulse-height distributions [7|]. 
A.G. Wright points out that the combination of the Poissonian nature of most 
light sources and the Binomial character of the photoelectric effect will give 
rise to a Poisson distribution for the number of detected photons. He goes on 
by establishing a method that uses this fact in an efficient way to build up the 
full pulse-height spectrum by using the response function of the corresponding 
detector. This method faces the fact that in many cases the response to a single 
photon is far from being Gaussian. Experimental input is used to accurately 
describe the amplitude distribution with no ad hoc assumption about how the 
detector will react to a single photon. 

In our MPHD model the violation of Poissonian statistics arises from the 
well known SiPM characteristics, incorporating the multiple-peak structure. 
It also makes use of the simplifying assumption that the individual peaks in 
the response function are Gaussian distributed. It was observed that when 
bias voltages are low and pixels are assumed to behave in a practically in- 
dependent way, the Gaussian hypothesis reproduces with great accuracy the 
measured spectra. Fig. 2(a)| shows the measured response function of a Pho- 
tonique devic^ZI for a bias voltage of 17.5 V after irradiation with a fluence 
of 3.1 x 10 9 electrons of 14MeV energy. Fig. 2(b)| shows the result of a ran- 
dom integration of noise signals for the same gate width. The relatively large 
number of counts outside the pedestal peak is a result of the high rate, r, of 
dark pulses. A simple convolution would add those events several times. SiPM 
are evolving toward the situation in which this algorithm can be successfully 
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Fig. 3. Independent clusters in which three pixels have been fired (crossed cells) 
in addition to the initial one (solid cell) as a consequence of optical cross-talk. 
Any of the possible 76 patterns can be derived from (a) to (1) by symmetry oper- 
ations. White cells represent pixels which remain off. The sum of all terms equals 
P(3) = <? 3 [36(1 - qf + 32(1 - qf + 8(1 - qf% where q is the probability of activa- 
tion of a single pixel. No distinction is made when there is more than one neighbour 
active. 

applied. Nowadays, devices are unfortunately not perfect and the high noise 
rate, for instance, would make the method inoperative. 

Optical cross-talk is modeled by considering the probability, g, of simultaneous 
activation of two isolated pixels when one of them has been triggered and 
deducing the probability for all possible clusters in the SiPM pixel matrix. 
No distinction is made when there is more than one neighbour active. Fig. [3] 
shows all independent cluster types for the case of three additional pixels. Only 
pixels sharing one side are allowed as part of a cluster. Cluster probabilities are 
given by the zero pixel cross-talk probability P(0) = (1 — g) 4 , and the TV-pixel 
cross-talk probabilities P(l) = 4g(l - qf, P(2) = g 2 [6(l - qf + 12(1 - qf], 
and P(3) = g 3 [32(l - qf + 32(1 - qf + 8(1 - q) 10 ]. The terms with (1 - q) 
are a consequence of the requirement that the pixels outside the cluster have 
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Fig. 4. Signals generated by the MPHD model. The straight line is the variable 
threshold level. Output of the leading edge discriminator is shown below. Superpo- 
sition of signals and optical cross-talk are visible in the first two peaks, a variation 
in the amplitude is observed for the next three single-pixel signals an after-pulse is 
visible in the next group of piled-up signals. 

to remain off. In absence of optical cross-talk, noise and after-pulses the pixel 
number distribution should be Poissonian. The mean value of this distribution, 
A, gives a measure of the photon detection efficiency. 

A realistic leading edge discriminator was implemented taking into account 
the blocking time of the module while the output is active (10 ns). Fig.[4]shows 
some of the pulses generated by the MPHD model, the threshold line and the 
output of the leading edge discriminator. Superposition of signals and optical 
cross-talk are visible in the first two peaks, a variation in amplitude can be 
observed in the next three single pixel signals and an after-pulse is visible in 
the next group of piled-up signals. 

Fig. [5] shows the measured pulse-height spectra for low amplitude signals in- 
duced by a short pulsed UV laseiliJ exciting a plastic scintillator. The data 
was taken with five different bias voltages. Underlying curves are the results 
of the MPHD model. The residuum between MPHD model and data is shown 
in panel (f), where the curves for difference bias voltages are off-set by nx 
0.004 to improve visibility. Table [I] summarises the resulting set of parame- 
ters. Fig. [6] shows the measured noise rate for a bias voltage of 17.9 V. The 
step structure is due to the multi-pixel events caused by optical cross-talk. 
Its sharpness is governed by the gain uniformity over the pixels and by the 
narrowness of the single pixel response function. Last steps are less defined 
due to the higher number of pixels involved. The pixels are better resolved 
with increasing: bias voltas es as a result of well known increase in the charge 
gain. The strong increase of noise rate with bias voltage is also visible. The 
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Fig. 5. Measured pulse-height spectra for five different bias voltages. Underlying 
curves are the result of the MPHD model described in the text, see Table [T] for the 
resulting set of parameters. Panel (f) shows the residuum between model and data, 
where the curves for difference bias voltages are off-set to improve visibility. 



open symbols are a result of the Monte Carlo. Simulated noise amplitude and 
gain variation have the same value as given in Table [U 
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Table 1 

MPHD model results for five bias voltages. The free parameters of the simula- 
tion are: mean signal amplitude (A), gain variation (cr^), dark pulse rate (r), op- 
tical cross-talk probability (p op t), after-pulse probability (p a ft), trap lifetime (r), 
mean number of detected photons (A), pedestal position (x pe d), and noise ampli- 
tude (a ped ). 
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Fig. 6. Measured noise rate as a function of discriminator threshold for a bias voltage 
of 17.9 V. The open symbols are the result of the MPHD model. Simulated noise 
amplitude and gain variation have the same values as given in Table [TJ 

3 SiPM irradiation 



A large increase of dark rate has been reported for relatively small radiation 
doses 1 81. After irradiation the noise rate increases because of the introduc- 
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tion of generation centers. This is the same behavior as for the bulk leakage 
current in diodes. In addition the creation of trapping centers is expected to 
increase the after-pulse probability. The application of the MPHD model de- 
scribed in the last section for the extraction of the characteristic parameters 
is considerably simpler when the multi-photoelectron peak structure is easily 
distinguishable. On the other hand the complete loss of this peak structure 
defines the limit of the application of SiPM for low light levels detection. The 
first irradiation dose was chosen so that only a small distortion of the peak 
structure was obtained. 14MeV electrons were used to irradiate a sample 
of SSPM-0701BG-TO18 diodes. The beam current was 10 nA. The electrons 
crossed a 0.3 mm thick aluminum window at 15 cm distance from the 1mm 2 
active area of the SiPM. Fluences on the diodes ranged from 3.1 x 10 9 to 
3.8 x 10 10 electrons/mm 2 . Heat dissipation and damage on the transparent 
epoxy layer protecting the silicon material were calculated and proved to be 
negligible. Grounding of the diode was provided in order to avoid damage by 
the sudden release of accumulated charge. 



4 Characterisation of the irradiation damage 

The SiPM were characterised before and after irradiation by studying the noise 
rates and the pulse-height spectra for low amplitude signals. Fig. [3 shows the 
noise rate before and after irradiation with 31 x 10 8 electrons of 14 MeV energy 
as a function of threshold in a leading edge discriminator. The figure includes 
curves for three different bias voltages. Two observations were made after the 
irradiation: 

(1) the rates of dark pulses are significantly larger. 

(2) the steps are much less pronounced than before irradiation. 

The simulation shows that an increase in the dark count rate is insufficient to 
explain the curves and it confirms that either an increase in the noise ampli- 
tude or the loss of gain uniformity (or a combination of both) can reproduce 
the measured values. 

The histograms in Fig. [8] show measured ADC spectra for bias voltages of 17.0 
and 17.4 V (compared to 17.9 V recommended voltage). The MPHD model was 
tested on this low photon yield data. The simulated integration of signals was 
performed in an interval of 30 ns, and the mean width of the single-pixel sig- 
nal was 9 ns. The pedestal peak in this case was clearly separated from the 
one-pixel peak. After-pulse probability was set to zero and the values for the 
noise rate and optical cross-talk probability were deduced from the measured 
noise rates as a function of discriminator threshold. The value of q was given 
to a good approximation by the ratio of noise rates above and below the first 
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Fig. 7. Measured noise rate before and after irradiation with 3.1 xlO 8 electrons/mm 2 
of 14MeV energy as a function of threshold in a leading edge discriminator. The 
plot shows curves for three different bias voltages. The step structure is due to the 
multi-pixel events caused by optical cross-talk. After irradiation the steps are less 
pronounced because of the much higher noise rate. 
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Fig. 8. Pulse-height spectra after irradiation with 3.1xl0 8 electrons/mm 2 of 14MeV 
energy for low amplitude signals. The underlying curves are a result of the MPHD 
model. 

step P(2)/P(l) = 4(g 3 — 2q 2 + q). The dark count rate might naively be taken 
from the measurement at lowest discriminator threshold, however, this value 
is not reliable due to the high probability of signal pile-up. If one assumes 
that the optical cross-talk probability was not significantly modified by the 
low radiation doses the dark count rate can be calculated by multiplying the 
noise rate for a threshold between one and two pixels by 1/4 (g 3 — 2q 2 + q). The 
distributions were well reproduced when noise amplitude and gain variation 
were allowed to change. The values for the gain variation before and after irra- 
diation were 0.08 chn and 0.12 chn, respectively. The noise amplitude changed 
from 1.49 chn before to 2.44 chn after irradiation. The good matching of the 
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Fig. 9. (a) Pulse-height spectrum obtained with larger laser intensity. Individual 
peaks are no longer visible. The ADC channel corresponding to the maximum proba- 
bility is a good measure of photon detection efficiency, (b) Relative photon detection 
efficiency for increasing electron fluence. 

two curves is a confirmation of the hypothesis of low after-pulse probability. It 
seems necessary to conclude that there has been an increase in leakage current 
and a severe loss of gain uniformity. 

For higher fluences the noise increase is so large that a multi-photoelectron 
peak differentiation is no longer possible. A different approach was used for 
dealing with such a situation: noise should not be relevant for large signals but 
a variation in the gain or photon detection efficiency should still be notice- 
able by studying the position of the maximum in the pulse-height spectrum, 
see Fig. [9](a). The result of such an investigation is shown in Fig. [9](b). The 
single pixel signal was monitored after each irradiation and no change in its 
amplitude was observable. The conclusion has to be that there has been a 
progressive reduction of the photon detection efficiency. This fact can be at- 
tributed to the loss of a progressively larger amount of pixels which remain 
permanently in the off state or to the reduction in the photon detection ef- 
ficiency of each pixel. Any of these two effects would be equally problematic 
for low light level detection. 



5 Conclusions 

Multi-anode photomultipliers are routinely used for scintillating fibre detector 
read-out showing all the drawbacks of conventional phototubes: stiff cables, 
high voltage power supplies and magnetic field sensitivity in addition to the 
specific problem of optical cross-talk among their many channels. If SiPM 
could be used instead an overall price reduction and a considerable detector 
simplification would be obtained. Conventional experimental techniques used 
to deal with high noise rates such as coincidences of several detectors can be 
complemented with more sophisticated methods based on intelligent trigger 
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algorithms implemented in FPGA chips in order to obtain a reliable tracking 
detector. Cooling Peltier modules would reduce dramatically the noise rates 
in case the mentioned methods prove to be insufficient [9[. 

For the Kaos/AI spectrometer at the Mainz Microtron, Germany, the viability 
of a fiber detector is studied as part of the electron arm tracking system in 
which SiPM are considered as possible photon detectors [10]. Simultaneously, a 
scintillating fibre tracker with SiPM read-out is being considered for the time- 
of-flight start detector or for a secondary active target in the future PANDA 
experiment [11] at FAIR. Here, radiation hardness is an important issue due 
to the relatively small distance from the target area to the detector position. 

A fluence of only 1.7 x 10 10 particles reduces already the number of detected 
photons by a factor of two. Good shielding will be necessary in many applica- 
tions where these doses are accumulated in a relatively short time. 
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